Genome Biology and Evolution

Animal mitochondrial rRNAs are commonly viewed as structurally reduced, yet sponge mt SSU rRNAs range from compact to highly expanded structures. Using nine conserved structural anchors, we compared 216 taxonomically resolved records from four classes and 22 orders, including 16 freshwater Spongillida and 200 marine sponges. Twelve homologous hypervariable substructures were coded as structural types, and their ordered combinations as composite types. We identified 38 structural types and 62 composite types across molecules ranging from 853 to 2,019 nt. Hexactinellida and freshwater Spongillida were each uniform for a distinct composite type but differed markedly in overall structure: hexactinellid mt SSU rRNAs were compact, whereas those of Spongillida were long and contained four to five candidate insertion regions. These results show that a conserved scaffold can accommodate extensive lineage-associated structural variation and provide a practical framework for comparing highly divergent mitochondrial rRNAs.

Chromosome evolution varies widely across Aphidomorpha: several aphid lineages exhibit extensive interautosomal reshuffling, whereas comparisons involving grape phylloxera and Adelges suggest greater chromosome conservation outside Aphididae. However, all published adelgid genomes represent the genus Adelges, leaving conservation across deeper adelgid divergences unresolved. Here, we report a chromosome-scale genome for the pine bark adelgid, Pineus strobi, providing the first genome for the genus Pineus and extending genomic sampling to an early-diverging adelgid lineage. Synteny analyses revealed broad conservation of major linkage groups between Pineus and Adelges across approximately 90 million years. The five grape phylloxera chromosomes also showed broad correspondence to the ten adelgid chromosomes, consistent with a limited number of chromosomal fusions or fissions and relatively little exchange among major linkage groups. These findings strengthen evidence that such extensive interautosomal reshuffling is not characteristic of Aphidomorpha as a whole. We also recovered a complete 25.3 kb mitochondrial genome with expanded, repeat-rich noncoding regions and complete circular genomes for the obligate nutritional symbionts "Candidatus Annandia pinicola" and "Candidatus Hartigia pinicola." Comparisons with Pineus similis symbionts revealed conserved coding capacity and genome-wide synteny, supporting conservation of this dual nutritional symbiosis in pine-associated adelgids. Finally, we recovered the first complete Wolbachia genome reported from an adelgid. Because Wolbachia can induce parthenogenesis, its presence in P. strobi raises the possibility that it contributes to the maintenance or reinforcement of parthenogenesis in a species lacking a viable sexual generation. Together, these genomes provide an integrated resource for aphidomorph evolution and symbiosis. Significance statementSeveral aphid lineages exhibit extensive reshuffling among autosomes, whereas previous comparisons with grape phylloxera and Adelges cooleyi suggested greater chromosome conservation outside Aphididae; however, all published adelgid genomes have represented Adelges, leaving the long-term stability of chromosome organization across Adelgidae unresolved. The first chromosome-scale genome from Pineus reveals broad conservation of major linkage groups across approximately 90 million years of adelgid evolution and clear correspondence with the five chromosomes of the grape phylloxera, supporting chromosome conservation outside Aphididae. The accompanying mitochondrial and symbiont genomes provide an integrated resource for studying aphidomorph genome evolution and symbiosis.

Movement of genetic material between non-parental organisms, called horizontal gene transfer (HGT), is well recognized in prokaryotes but represents an underestimated force for the acquisition of novel traits in eukaryotes. The mechanisms of cross-domain HGT remain poorly understood despite numerous reports of its occurrence. Thus, progress in the field remains limited by the lack of targeted approaches for detecting HGT events. Using BLASTp search and sequence identity, we describe an HGT-derived protein from the chytridiomycete fungus Neocallimastix californiae. This protein has two cysteine-rich repeats (CysRReps) and other functional domains that are highly similar to those of prokaryotes. Based on the alignment of several CysRReps, we identified 859 additional eukaryotic proteins spanning 43 protein families, each with at least one match to bacterial, archaeal, or viral proteins. HGT-derived proteins with CysRReps belong to multiple eukaryotic taxa, some of which are well-known HGT models. Bacterial, archaeal, and viral proteins comprise the prokaryotic counterparts. Bacteria belong to 15 phyla, spanning a diverse array of physiologies, habitats, and lifestyles. Viral diversity is largely restricted to Caudoviricetes, suggesting a virus-mediated mechanism of DNA integration underlying these transfers. Although the function of CysRReps remains unclear, their association with HGT events across diverse eukaryotic and prokaryotic taxa suggests the existence of a universal molecular mechanism that facilitates gene transfer across phylogenetic boundaries.

The phylogenetic position of Nanohaloarchaea has been debated, these nanosized archaea being alternatively proposed as sister group to Haloarchaea, members of the DPANN-Archaea, or sister group to Methanocellales. Screening a set of universal proteins, we identified four insertions located at critical locations in three ribosomal proteins and one RNA polymerase subunit that support the branching of Nanohaloarchaea as sister group to Aenigmarchaea within DPANN cluster II (sensu Dombrowski et al., 2020). Insertion analyses and phylogeny of the monomeric primase specific to DPANN-Archaea confirm the existence of a robust clade grouping Undinarchaea, Naiadarchaea and DPANN cluster II, that we propose to call Nanostetteria. Our insertion analysis also supports including Altiarchaea within DPANN-Archaea and suggest a new clades that has not been recovered in phylogenetic analyses, one grouping DPANN-Archaea with Stygia (Hadarchaea and relative) and an even large one grouping these lineages with Acherontia (Thermococci and relatives). The insertion defining this larger clade, present in the ribosomal protein uS7, is also present at the same position in Thaumarchaea, Korarchaea and a subgroup of Asgardarchaea. Whereas the insertion in Thaumarchaea is certainly due to an independent event, we discuss alternative hypotheses that can explain those present in Korarchaea and Asgardarchaea. Finally, we noticed several cases of MAGs misannotations, indicating that insertion analysis can be useful to identify protein with misleading affiliations. The existence of insertions in otherwise highly conserved universal proteins involved in translation or transcription could partly explain the high rate of protein evolution in some archaeal lineage, especially in DPANN-Archaea.

Genome sizes vary across four orders of magnitude in flowering plants, with consequences for evolution. Much of this variation is due to differences in transposable element content, particularly in long terminal repeat (LTR) retrotransposons. Despite their importance in plant genome evolution, LTRs have long been challenging to characterize because of their repetitive nature, but recent advances in sequencing allow more detailed explorations of their behavior. They are now known to occupy distinct genomic niches, and to evolve and proliferate over time as they escape host controls. In this study, we present a new genome sequence of the largest Brassicaceae genome, Hesperis matronalis, and describe its transposable element complement and how LTRs contributed to its genome expansion. We find evidence for both early proliferation of Ty3 elements and rapid recent expansion of Ty1-copia elements. In addition, we place the H. matronalis LTRs in a broader context of LTR evolution in the Brassicaceae, showing that the dominant copia families are part of an evolutionary radiation endemic to Hesperis. Finally, we describe differences in LTR age, proximity to genes, and apparent removal rate which suggest consistent genomic niches over the lifespan of a TE family. These results shed light on how LTRs evolve dynamically with host genomes, and have contributed to the expansion of the largest genome in the Brassicaceae. Significance StatementOrganisms differ in many ways, including the size of their genomes. Genome size, in turn, can affect many aspects of evolution, such as what kinds of mutations are available to natural selection and how effectively natural selection can act. The "complexity" of an organism does not predict its genome size; rather, much of this variation is explained by the amount of transposable elements, or "jumping genes," that inhabit a genome. This study describes the evolutionary history of the transposable elements in dames rocket (Hesperis matronalis), which has the largest genome of any plant in the mustard family (Brassicaceae). We find that the genome inflation in Hesperis is due to transposable elements that are not found elsewhere in the Brassicaceae, and likely diversified within Hesperis. Characterizing the evolution and behavior of transposable elements in this genome offers a better understanding of the forces that determine genome size and ultimately affect the evolution of life on earth.

Although numerous studies have documented horizontal gene transfer (HGT) from bacteria into insect genomes, they have been heavily biased toward the endosymbiont Wolbachia. In contrast, comparatively few studies have examined HGT from other symbionts and insect-associated bacteria. Moreover, research on these systems has primarily focused on transferred genes that retain function in the insect host, contributing to metabolism, symbiosis, detoxification, or host adaptation. Consequently, relatively few empirical systems have been available to investigate the persistence, degradation, diversification, and vertical inheritance of apparently non-functional HGT-derived sequences over macroevolutionary timescales. The widespread occurrence of putatively Arsenophonus-derived sequences across Cimicidae therefore provides a rare opportunity to extend these investigations beyond the Wolbachia model and to examine the long-term evolutionary fate of non-functional HGT-derived DNA across tens of millions of years of host diversification. We show that Arsenophonus-derived HGTs have multiple independent origins across Cimicidae but that one major HGT lineage has persisted through diversification of the subfamily Cimicinae for at least 50 million years. Phylogenetic and compositional analyses indicate that an ancestral Arsenophonus-derived genomic region has undergone progressive fragmentation, leaving numerous dispersed, apparently non-functional remnants while preserving a clear evolutionary signature. These results extend the study of bacterial HGT beyond the Wolbachia model and demonstrate that non-functional symbiont-derived DNA can persist over macroevolutionary timescales.

Whole-genome duplication (WGD; or polyploidy) has played a major role in the evolution of many lineages however, our understanding of the processes that shape genome evolution following WGD remains incomplete. While polyploidy duplicates the entire genome sequence, it is rediploidisation that establishes independent duplicated genes. Rediploidisation proceeds through suppression of meiotic recombination across polysomic loci thus restoring disomic inheritance, a process that is not synchronised across the genome. Despite its importance, the mechanisms underlying this process remain poorly understood. The slowly evolving genomes of paleopolyploid Acipenseriformes paddlefish and sturgeon provide an invaluable system for investigating this, as rediploidisation was highly asynchronous in these lineages. In both genomes ohnologs tend to segregate into blocks on the chromosomes according to rediploidisation timing, a pattern that suggests links between chromosomal structure and rediploidisation. Here, we analyse a newly-produced duplicate-resolved paddlefish genome assembly and show a strong concordance between genome rearrangement and rediploidisation timing. We also find that topologically associated domain (TAD) boundaries are associated with rediploidisation block boundaries. Together these results indicate that rediploidisation in acipenseriformes occurred through a process of genome rearrangements that was subject to functional constraints imposed by 3D genome architecture. We investigate the evolution of Hox clusters in these lineages, revealing a previously overlooked duplicate HoxC region in paddlefish, and both ancestral and lineage-specific Hox cluster rediploidisation with substantially different timings. These findings highlight a complex evolutionary history following WGD in Acipenseriformes with implications for understanding short-term adaptations to polyploidy as well as longer-term diversification of lineages.

Background: Chromosome-level genome assemblies are increasingly enabling tests of chromosome evolution, conserved synteny, and sex chromosome conservation across diverse animal lineages, including spiders. Results: Here, we present a chromosome-level genome assembly for the African hermit spider, Nephilingis cruentata, a species with extreme female-biased sexual size dimorphism and a cytogenetically inferred XX2 sex chromosome system. The final Hi-C-assisted assembly spans 1.72 Gbp, with 99.5% of bases assigned to 13 pseudochromosomes, a scaffold N50 of 131.6 Mbp, and a BUSCO completeness score of 98.8%. We annotated 20,021 protein-coding genes, and repetitive elements accounted for 42.7% of the genome. Sex-specific whole-genome resequencing identified Chr02 and Chr07 as candidate X chromosomes based on reduced male coverage, consistent with the expected XX2 system. Using comparative whole-genome alignments across existing chromosome-scale spider assemblies, we also show that sex-linked chromosomes retain broad homologous identity across sampled spider lineages but exhibit lower synteny conservation and greater chromosome-length divergence than autosomes. Conclusions: These results suggest that spider sex chromosomes are conserved in homologous identity but more labile in structure, providing a comparative framework for studying sex chromosome conservation and divergence across Araneae.

The eriophyoid mites (Acari: Eriophyoidea) represent an extreme case of genome streamlining, with genomes averaging just 32Mb, among the smallest of all animals. This clade of mites is highly diverse with more than 4000 species. Their morphologies are specialized for feeding on their host plants, including a simplified worm-like body plan with just two pairs of legs and modified mouth parts that can induce the formation of galls or plant deformities during feeding. Their generally strong host affinities make eriophyoid mites appealing for use as biological control agents, although their short generation times and small genomes could facilitate rapid evolution and impact their efficacy in management programs. Here, we sequenced the genome of the biological control mite Floracarus perrepae, producing a highly contiguous genome totalling just 23.1 Mb, among the smallest of all animals. We also assembled another non-eriophyid mite genome from accidental DNA bycatch (69.4Mb). We placed this new genomic resource in a phylogenetic context to reveal that mites have highly dynamic genome evolution, with a significant trend in genome downsizing in the eriophyoids. Our results suggest that this streamlining is associated with non-genic elements such as the suppression or excision of retrotransposons and purging of introns. As new sequencing techniques become available, novel genomic resources for tiny organisms such as F. perrepae will be more readily accessible, facilitating both fundamental genome evolutionary biology and applied sciences such as biological control programs which use eriophyoid mites for the management of invasive species.

Transposable elements (TE) jump from one genomic locus to another. Since increase in their copy number is a metabolic burden for the host, TE are considered as genomic parasites. Although host-TE co-existence is regarded as an evolutionary arms race, the hypothesis is not extensively tested especially using evolutionary genomics. We provide a hypothesis testing framework to understand the distribution of TE in genic regions of the host genome, variation in the regulation of TE by host, and effect of these two factors on host-TE co-evolutionary dynamics. We test our hypothesis by understanding the distributions of potentially active TEs in the genome of 78 teleost fishes, representing major families and orders within the clade. Our analysis reveals coevolutionary arms race predicted by the Red Queen dynamics.

The European leaf-toed gecko (Euleptes europaea) is a small, nocturnal gecko endemic to the western Mediterranean. As a phylogenetically distinctive member of the Gondwanan family Sphaerodactylidae, it represents an important species for studying Mediterranean island biogeography, adaptation, and reptile genome evolution. The species also occupies a key position for investigating the evolution of sex chromosomes, as geckos exhibit remarkable diversity and frequent transitions in sex-determination systems. We present a chromosome-level genome assembly of Euleptes europaea generated as part of the Vertebrate Genomes Project. The 1.8 Gb assembly has a scaffold N50 of 102.3 Mb (contig N50 27 Mb), with 21 chromosome-scale scaffolds corresponding to the known karyotype (2n = 42). The primary assembly has a BUSCO completeness of 97.80% (95.60% as single-copy), a k-mer completeness of 96.00%, and a k-mer quality value (QV) of 61.20. Repetitive elements account for 53.20% of the genome and genome annotation identified 18,633 protein-coding genes. This high-quality reference genome will facilitate studies of genome evolution, island adaptation, and sex chromosome evolution across geckos and other reptiles.

Most cyanobacterial genomes are rich in the GCGATCGC octamer, also known as Highly Iterated Palindrome 1 (HIP1). Despite its description over three decades ago, the biological function of this highly abundant sequence is only beginning to be elucidated. HIP1 is recognized by two DNA methylases, DmtA and DmtC, and is characterized by its evolutionary conservation and a quasi-periodic distribution within genomes. However, whether the phylogenetic distribution of HIP1 correlates with the presence of functional categories of protein families remains unknown. Here we investigated whether certain protein families share a phylogenetic distribution with this abundant palindromic sequence across cyanobacterial genomes. Our analysis revealed a strong phylogenetic correlation between several proteins of the Type IV secretion system (T4SS) and the abundance of HIP1. This finding aligns with recent discoveries demonstrating that HIP1 enhances DNA transformation in a methylation-dependent manner in two distinct cyanobacterial species. Consequently, we hypothesize that HIP1 function as a conserved adaptation for horizontal gene transfer (HGT) at the phylum level, potentially by serving as a DNA-uptake recognition sequence in cyanobacteria. Significance statementScientists have long been baffled by the HIP1 sequence, a short, highly common, repetitive DNA pattern found across almost all cyanobacterial genomes. Our study used a whole-genome evolutionary approach and found that the presence of this repetitive pattern is tightly linked to the presence of a cells external DNA uptake system. This tight co-evolutionary relationship suggests that HIP1 isnt just random genomic feature, but a conserved evolutionary adaptation used by the entire cyanobacterial phylum to specifically enhance their ability to acquire new genes from one another.

The coffee bean weevil, Araecerus fasciculatus (Coleoptera, Curculionoidea, Anthribidae), is a cosmopolitan pest of over 100 stored agricultural commodities, with particular economic impact on coffee (Coffea arabica). Although two chromosome-level anthribid genomes have recently been released as part of the Darwin Tree of Life (DToL) project (Booth et al. 2024; Crowley et al. 2025), no functionally annotated genome has been available for the family. Here we present a draft genome assembly for A. fasciculatus, generated from PacBio HiFi long reads and processed through a three tiered metagenomic filtering pipeline to remove host plant (C. arabica) and microbial contamination. The final assembly spans 475 Mb across 3,617 scaffolds (N50 = 170 kb) with 88.5% BUSCO completeness (insecta_odb10) and only 3.1% duplication. Gene prediction with BRAKER2 identified 22,384 protein-coding genes, of which 11,783 received functional annotations through SwissProt similarity. Notably, we identified 92 cytochrome P450 (CYP) genes, including tandem gene clusters on two scaffolds (4 genes on ptg000464l, 5 genes on ptg001867l), suggestive of lineage-specific expansion through tandem duplication. Homology searches against Drosophila melanogaster caffeine-metabolizing P450s (CYP12D1, CYP6d5, CYP6a8) recovered strong matches (e-values 9.7 x 10-110 to 5.4 x 10-101, 33-38% identity). In stark contrast, comprehensive BLAST searches for bacterial caffeine N-demethylase genes (ndmA/B/C/D), which mediate caffeine degradation via horizontal gene transfer in the coffee berry borer Hypothenemus hampei (Scolytinae), returned zero hits across the A. fasciculatus genome, predicted proteome, and associated bacterial scaffolds. AlphaFold2 structure prediction of four top Araecerus P450 candidates produced high-confidence models (pLDDT 84.5-93.9, pTM 0.735-0.930) with conserved P450 catalytic motifs. Foldseek structural homology searches confirmed that all four candidates adopt cytochrome P450 folds (top hits: human CYP3A4, CYP3A7, CYP11A1; TM-scores 0.90-0.92; probability 1.000), with zero hits to bacterial Rieske-fold enzymes. Molecular docking of caffeine against these structures yielded binding affinities of -5.41 to -5.80 kcal/mol for the Araecerus candidates, comparable to or exceeding the -5.55 kcal/mol obtained for the experimentally validated Drosophila CYP6a8 and substantially stronger than the -3.70 kcal/mol for the bacterial NdmA structural outgroup (PDB: 6ICP). Phylogenetic analysis revealed that all four candidates have clear orthologs in two non-seed-feeding DToL anthribids (Pseudeuparius sepicola and Platystomos albinus), demonstrating that these P450 genes predate the dietary transition to caffeine-containing seeds. The Araecerus candidates predominantly belong to the CYP6 family (clan 3), whereas the primary Drosophila caffeine P450 CYP12D1 belongs to the mitochondrial clan, confirming convergent recruitment of different P450 subfamilies for caffeine metabolism. These results support the hypothesis that A. fasciculatus employs an insect-encoded, P450-mediated caffeine detoxification pathway fundamentally distinct from the bacterial horizontal gene transfer mechanism documented in Scolytinae. This represents convergent evolution of caffeine resistance via independent molecular strategies within Curculionoidea, and provides the first functionally annotated genomic resource for comparative studies across the Anthribidae.

The rat tapeworm, Hymenolepis diminuta, is an important laboratory model for uncovering molecular processes that underly the success of tapeworms as parasites. Despite its importance, a high-quality reference genome for this species is lacking. Here we present a highly contiguous and effectively complete genome of H. diminuta assembled from PacBio HiFi long-read sequencing data. Our primary assembly consists of 7 scaffolds (N50=29.25 Mb) with total length of 186.53 Mb, has only 7 gaps, and contains 95.7% complete Lophotrochozoan BUSCOs. Our assembly allows us to confirm aspects of Hymenolepis genome organization, such as high repeat content and unusual chromosomal ends, and to show that Hymenolepis genomes encode [~]10,000 genes. Together with annotations of nuclear tRNAs, mtDNA protein coding genes, and mtDNA tRNAs, our assembly currently provides one of the most complete genome resources for a tapeworm species and will enable research on parasitism, animal regeneration, development, and evolution.

Rates of de novo mutations ({micro}) are highly variable among different taxa but less is known about their variability at the intraspecific level. Using a large data set (n=364 trios) from eight nine-spined stickleback (Pungitius pungitius) populations differing in their effective population sizes (Ne), we tested the prediction of the drift-barrier hypothesis (DBH) that {micro} scales negatively with increasing Ne. Indeed, {micro} was a negative function of Ne, also after correcting for phylogenetic non-independence of populations. While the range of variation in mutation rates across populations spanned more than a 4-fold range ({micro} = 1.60 - 6.99 x 10-9), we discovered one highly mutable family with a five-times higher mutation rate than the population average. Evidence was also found (i) for reduced efficiency of selection in small freshwater populations subject to strong genetic drift, (ii) that maternal age associates positively with {micro}, and shorter generation time elevates per-year mutation rates, (iii) that both replication errors and DNA repair efficiency contributed to {micro}, and that (iv) mutation rate variation has a polygenic basis. In general, the results provide support for the DBH identifying elevated mutation rates in populations subject to strong genetic drift.

Gene duplication is a major source of genetic variation and is considered as an important driver of evolutionary novelties, including eusociality in insects, such as ants, wasps, bees, and termites. However, it remains unclear whether selection acts to increase gene copy number in social insects. Here, we studied the paranomes, namely, the whole set of paralogous genes in a genome, of Blattodea and Hymenoptera. We estimated the rates of gene duplication and loss using the distribution of synonymous substitution rate of paralogous genes and showed that regardless of sociality, duplicated genes were lost more rapidly than expected under random drift, indicating that negative selection on duplicated genes is prevalent across Blattodea and Hymenoptera. The rates of gene duplication and loss varied independently of sociality levels, but were positively related to genome size, suggesting the expansion/contraction of gene families can be a side effect of genome expansion/contraction. These results call for a reevaluation of adaptative gene duplications.

Throughout eukaryotic evolution, the structure of the ribosome has been highly conserved, featuring 80 common protein gene families. However, in many eukaryotes, paralogs of these proteins are present. "Specialised ribosomes" have been documented across diverse groups of eukaryotes where they play an important role in the regulation of translation of specific mRNAs. In the case of specialised ribosomes it has been documented that assembled ribosomes that contain specific paralogs can directly affect translational output. This has been proposed to contribute to the regulation of complex responses to environmental change and to coordinate cell-type specific physiology. This poses the question of whether ribosome specialisation principally emerges under an adaptive or neutral model of evolution. Using gene tree-species tree reconciliation, we test competing hypotheses regarding the evolutionary drivers of ribosome specialisation. We determine that examples of specialisation tend to emerge by independent duplication of the same ribosomal proteins in different lineages. We show that pathways to specialisation through paralog formation have arisen independent of: (i) paralog location within the 3D ribosome complex, and (ii) positive selection in these paralogs. We determine that the generalisable model of best fit for the evolution of paralog-mediated eukaryotic ribosomal specialisation is one of constructive neutral evolution. In lineages with small effective population sizes and increased complexity, the emergence and retention of ribosomal protein paralogs has provided the raw material for ratcheting and the emergence of translational regulation at the level of the ribosome.

Multicellularity has evolved multiple times across the eukaryotic tree of life, including among protist lineages. Because transposable elements (TEs) strongly influence genome architecture and gene regulation, understanding their potential impact on genome structure and their relationship with gene expression may provide insight into the evolution of multicellularity. Here, we generated a new genome assembly for the facultatively multicellular amoeba Acrasis kona and performed comparative analyses of TE composition, TE diversity, and TE-density organization across diverse protist lineages. Comparative analyses included unicellular and multicellular representatives from across the tree of eukaryotes, (Heterolobosea, Filasterea, Cristidiscoidea, and Chlorophyceae), including Naegleria spp., Tetramitus jugosus, Capsaspora owczarzaki, Pigoraptor spp., Fonticula alba, Parvularia atlantis, Volvox carteri, and Chlamydomonas reinhardtii. To examine relationships between TEs and gene regulation, we integrated transcriptomic datasets from A. kona, Capsaspora owczarzaki, and Volvox carteri with genome-wide TE-density analyses of differentially expressed genes. TE abundance and composition varied substantially among lineages, with species that exhibit more complex developmental or cellular organization generally containing higher TE proportions than closely related unicellular taxa. Patterns of TE-density organization near up-regulated, down-regulated, and non-differentially expressed genes also differed among systems, ranging from strong TE depletion in A. kona to weaker or cell-type-specific patterns in Capsaspora and Volvox. Together, these findings suggest that transposable elements are associated with multicellularity across diverse protist lineages, although the specific roles they play appear to be complex, lineage-specific, and not yet fully understood.

The evolutionary transition to terrestrial life required overcoming several physiological hurdles; however, such challenges were amplified in desert environments. While several xeric-adapted arachnids utilize permanent burrows or "sit-and-wait" foraging strategies as possible energy conservation adaptations in harsh habitats, camel spiders exhibit a counterintuitive, high-energy lifestyle. To investigate the molecular underpinnings distinguishing Solifugae within Chelicerata, we utilized a comparative genomics framework that incorporates a newly sequenced, previously unpublished solifuge genome. We identified lineage-specific expanded orthogroups and evaluated selective pressures acting upon paralogous sequences within our ingroup solifuge species. Additionally, we also focused on fatty acid-associated proteins and heat shock proteins to elucidate how Solifugae may have evolved such anomalous behaviors compared to their arachnid relatives. Our analyses revealed significant signatures of positive selection within key gene families across the solifuge lineage. Notably, paralogs within the cytochrome P450 and biotinidase families showed consistent evidence of selection across all three taxa, suggesting specialized metabolic or detoxification requirements. Furthermore, we identified candidate loci implicated in axonal guidance and lipid metabolism, and a specialized fatty acid enzyme repertoire. While subsequent research is required to determine whether some of the genomic signatures unveiled here are shared across a broader phylogenetic distribution within Solifugae, we establish a critical baseline for future functional validation.

We present a high-quality chromosome-level reference genome for the Saimaa ringed seal (Pusa saimensis), an endangered freshwater pinniped endemic to Lake Saimaa, Finland. The assembly spans 2.353 Gb and comprises 15 autosomes together with the X and Y sex chromosomes. Using Oxford Nanopore Technologies (ONT) long-read sequencing and Hi-C scaffolding, we achieved a telomere-to-telomere assembly for all chromosomes, except the Y chromosome. Genome annotation identified approximately 21,800 protein-coding genes, consistent with other mammalian genomes. Assembly completeness was high, with BUSCO analysis recovering 99.6% of expected complete single-copy genes (98.2% single-copy and 1.3% duplicated). Comparative analyses revealed a highly conserved chromosomal architecture, with only minor syntenic differences relative to other pinniped chromosome-level assemblies. Previously described cytogenetic fusion events in Phocidae were confirmed (chromosomes 2 and 7). A translocation between chromosomes 6 and 7 distinguishes phocids from the otariids. In general, more distantly related taxa exhibit an increasing degree of intrachromosomal rearrangements. Notably, we identified a large intrachromosomal rearrangement on chromosome 2 that appears specific to the Saimaa ringed seal. Phylogenomic analysis based on 9,226 single-copy orthologues placed the Saimaa ringed seal as a sister lineage to the Baltic ringed seal (Pusa hispida botnica), while confirming also other established evolutionary relationships among pinnipeds. Comparative gene family analysis between the Saimaa ringed seal and the closely related grey seal (Halichoerus grypus) revealed lineage-specific differences driven by a limited number of gene families. In the Saimaa ringed seal, expansions were observed in ion transport, cytoskeleton, and regulatory genes, potentially reflecting adaptation to freshwater conditions. In contrast, the grey seal showed expansions in olfaction, immune-and spermatogenesis-associated gene families, including MAGE/MIA genes, consistent with differences in ecology and mating systems. This reference genome provides an important resource for studies of pinniped genome evolution, as well as conservation and population genomics of the Saimaa ringed seal, facilitating future work on genetic diversity, inbreeding, mutational load and adaptive potential in this highly endangered species.